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(Dated: September 19, 2006) 

We report on a measurement of the B d mixing frequency and the calibration of an opposite-side fla- 
vor tagger in the D0 experiment. Various properties associated with the b quark on the opposite side 
of the reconstructed B meson were combined using a likelihood-ratio method into a single variable 
with enhanced tagging power. Its performance was tested with data, using a large sample of recon- 
structed semileptonic B — ► fiD° X and B — * fiD* X decays, corresponding to an integrated luminos- 
ity of approximately 1 fb _1 . The events were divided into groups depending on the value of the com- 
bined tagging variable, and an independent analysis was performed in each group. Combining the re- 
sults of these analyses, the overall effective tagging power was found to be eV 2 = (2.48±0.21+ ; 06 )%. 
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The measured B% mixing frequency Ami = 0.506 ± 0.020 (stat) ± 0.016 (syst) ps 1 is in good 
agreement with the world average value. 



I. INTRODUCTION 



(iVtot) that are tagged (iV tag ): 



Particle-antiparticle mixing in the B° (-8°) system has 
been known for more than a decade now [lj and has 
been studied at the CERN LEP collider and subsequently 
at the Fermilab Tevatron collider during Run I. It is 
currently being measured at the B-factory experiments, 
Belle and BaBar, and the Fermilab Tevatron collider ex- 
periments during Run II. 

Mixing measurements involve identifying the "flavor" 
of the B° meson at production and again when it decays, 
where "flavor" indicates whether the meson contained a 
b or a b quark. The decay flavor is identified from the B° 
decay products when the B° meson is reconstructed. The 
determination of the initial "flavor" is known as flavor 
tagging. 

The B® meson flavor at its production can be identified 
using information from the reconstructed side or from 
the opposite side (see Fig. Q). One can tag the flavor 
using charge correlation between "fragmentation tracks" 
associated with the reconstructed B meson. Such corre- 
lations were first observed in e + e~ — > Z a — ► bb events by 
the OPAL experiment 0- This is known as "same-side 
flavor tagging." The flavor can also be inferred from the 
decay information of the second B meson in the event, 
assuming that b and b are produced in pairs, and thus in 
the ideal case, the two mesons have opposite flavors. This 
method is known as "opposite-side flavor tagging." An 
advantage of the latter method is that its performance 
should be independent of the type of the reconstructed 
B meson. 

Measurement of the B^ mixing parameter is an impor- 
tant test of the opposite-side flavor tagging as the same 
tagger is used for our study of B s mixing. Studies of 
tagged B° and B + samples at hadron colliders could re- 
veal physics beyond the standard model |j] . Finally, this 
technique of flavor tagging developed at the Tevatron can 
also be useful for future experiments at the Large Hadron 
Collider at CERN. 

This paper describes the opposite-side flavor tagging 
algorithm used by the D0 experiment in Run II and the 
measurement of its performance using B — > X and 

B — ► fi + D*~X events. Throughout the paper, a reference 
to a particular final state also implies its charge conju- 
gated state. B + decays represent the main contribution 
to the B — > n + D°X sample, and B° decays dominate 
in the B — > n + D*~X sample. We measure the flavor 
tagging purity independently for reconstructed B + and 
B° events and then extract the B° oscillation frequency. 
This technique allows us to verify the assumption of inde- 
pendence of the opposite-side flavor tagging on the type 
of reconstructed B meson. Its performance is described 
by the two parameters, efficiency and dilution. The effi- 
ciency e is defined as the fraction of reconstructed events 



JVtag/JVtot- 



(1) 



The dilution T> is a normalized difference of correctly 
and wrongly tagged events: 



V = 



Nr., 



N COT + Nv, 



2P-1, 



(2) 



lag 



where P = N COI /N tag is called the purity. The terms 
"correctly" and "wrongly" refer to the determination of 
the reconstructed B meson flavor. The effective tagging 
power of a tagging algorithm is given by eT> 2 . 



Reconstructed side 




PV 



B 



Opposite Side 



7C/K x 




FIG. 1: Diagram of an event with a reconstructed B° candi- 
date. PV indicates the primary vertex for the event. 



II. 



DETECTOR DESCRIPTION 



The D0 detector is described in detail elsewhere |||. 
The main features of the detector essential for this anal- 
ysis are summarized below. Tracks of charged particles 
are reconstructed from the hits in the central tracking 
system, which consists of the silicon microstrip tracker 
(SMT) and the central fiber tracker (CFT), both lo- 
cated in a 2 T superconducting solenoidal magnet. The 
SMT has 800,000 individual strips, with a typical 
pitch of 50 — 80 fim. and a design optimized for tracking 
and vertexing capability for \r)\ = 3. The pseudorapid- 
ity, r] = — In [tan(#/2)], approximates the true rapidity, 
V = \ \n[(E+p z c) / (E— p z c)}, for finite angles in the limit 
of (mc 2 /E) — ► 0, $ being the polar angle. We use the 
term "forward" to describe the regions at large | rj | . The 
SMT system consists of six barrels arranged longitudi- 
nally (each with a set of four layers of silicon detectors ar- 
ranged axially around the beam pipe), interspersed with 
16 radial disks. The CFT has eight thin coaxial barrels, 
each supporting two doublets of overlapping scintillating 
fibers of 0.835 mm in diameter, one doublet being parallel 
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to the beam axis and the other alternating by ±3° rel- 
ative to this axis. Light signals are transferred via clear 
light fibers to solid-state photon counters (VLPCs) that 
have ss 80% quantum efficiency. 

The muon system consists of a layer of tracking detec- 
tors and scintillation trigger counters in front of 1.8 T 
toroids, followed by two additional similar layers after 
the toroids. Muon tracking for |?/| < 1 relies on 10-cm- 
wide drift tubes, while 1-cm mini drift tubes are used for 
1<N<2. 

Electrons are identified using matching between the 
tracks identified in the central tracker and energy de- 
posits in a primarily liquid-argon/uranium sampling 
calorimeter We also use the energy deposits in 

the central preshower detector which consists of 
three concentric cylindrical layers of triangular scintilla- 
tor strips and is located in a nominal 5 cm gap between 
the solenoid and the central calorimeter, to provide ad- 
ditional discrimination between electrons and fakes. The 
calorimeter consists of the inner electromagnetic section 
followed by the fine and coarse hadronic sections. In this 
analysis, we only use the central calorimeter (| rj |< 1). 



III. DATA SAMPLE AND EVENT SELECTION 

This measurement is based on a large semileptonic 
B decay data sample corresponding to approximately 
1 fb _1 of integrated luminosity collected with the D0 
detector between April 2002 and October 2005. 

B mesons were selected using their semileptonic de- 
cays B — ► [i + vD°X and were divided into two exclu- 
sive groups: the D* sample, containing all events with 
reconstructed D*~ — > D°ir~ decays, and the D° sample, 
containing all the remaining events. The D* sample is 
dominated by B® — * [i + vD*~X decays, while the D° 
sample is dominated by B + — > /.i + vD° X decays. 

The flavor tagging procedure was developed using 
events from the D° sample. Events from the D* sample 
were used to measure the purity of the flavor tagging and 
the oscillation parameter Am^. In addition, the purity 
was measured in the D° sample to test the hypothesis 
that the flavor tagger is independent of the type of re- 
constructed B meson. 

Muons for this analysis were required to have hits in 
more than one muon chamber, an associated track in 
the central tracking system with hits in both SMT and 
CFT detectors, transverse momentum p T > 2 GeV/c, as 
measured in the central tracker, pseudorapidity < 2, 
and total momentum p M > 3 GeV/c. 

All charged particles in a given event were clustered 
into jets using the DURHAM clustering algorithm [|| 
with the cut-off parameter set to 15 GeV/c. Events with 
more than one identified muon in the same jet or with 
the reconstructed J/tp — ► decays were rejected. 

D° candidates were constructed from two tracks of 
opposite charge belonging to the same jet as the recon- 
structed muon. Both tracks were required to have trans- 



verse momentum px > 0.7 GeV/c and pseudorapidity 
\rj\ < 2. They were required to form a common D vertex 
with a fit x 2 < 9, number of degrees of freedom being 1. 
For each track, the projection ex (onto the axial plane, 
i.e. perpendicular to the beam direction) and projection 
6l (onto the stereo plane, i.e. parallel to the beam di- 
rection) of its impact parameter with respect to the pri- 
mary vertex, together with the corresponding uncertain- 
ties (cr(er), cr(ei)) were computed. The combined impact 

parameter significance S — \f[er / ,(7 (^t)] 2 + [(-l/^I^l)] 2 
was required to be greater than 2. The distance be- 
tween the primary and D vertices in the axial plane was 
required to exceed 4 standard deviations: d^ /cr(d^) > 4. 
The accuracy of the d® determination was required to 
be better than 500 fj,m. The angle a® between the D° 
momentum vector and the direction from the primary 
to the D vertex in the axial plane was required to sat- 
isfy the condition cosa^ > 0.9. The tracks of the muon 
and D° candidate were required to form a common B 
vertex with a fit x 2 < 9, with number of degrees of free- 
dom being 1. The mass of the kaon was assigned to the 
track having the same charge as the muon; the remain- 
ing track was assigned the mass of the pion. The mass 
of the (fi + D° ) system was required to fall within the 
2.3 < M(fi+D° ) < 5.2 GeV/c 2 range. 

If the distance g?|? between the primary and B vertices 
in the axial plane exceeded 4cr(d|?), the angle a|? between 
the B momentum and the direction from the primary to 
the B vertex in the axial plane was required to satisfy the 
condition cosa|? > 0.95. The distance d|? was allowed to 
be greater than d®, provided that the distance between 
the B and D vertices d^ D was less than "iaid^P). The 
uncertainty a(d^) was required to be less than 500 /xm. 
In addition, the cut pt(D° ) > 5 GeV/c 2 was applied. 

To select [i + D*~~ candidates, we searched for an addi- 
tional pion track with px > 0.18 GeV/c and the charge 
opposite to the charge of the muon. The mass differ- 
ence AM =_M(Ktttt) - M(Ktt) for D* candidates, with 
1.75 < M(D°) < 1.95 GeV/c 2 , is shown in Fig.|2 The 
peak corresponding to the mass of the soft pion in the 
H+D*- sample is clearly seen. 

All events with 0.1425 < AM < 0.1490 GeV/c 2 were 
included in the D* sample. The remaining events were 
assigned to the D° sample. The Kir mass distributions 
for these two samples together with the results of the fits 
are shown in Figs. |j3and|U The procedure to fit these 
mass spectra is described in Sec I VIII In total, 230551 ± 
1627 B -> fi+vD decays and 73532 ± 304 B -> fi+vD* 
decays were reconstructed. 



IV. VISIBLE PROPER DECAY LENGTH 

The oscillations of B mesons are usually studied as a 
function of their proper decay length. Since in semilep- 
tonic B decays an undetected neutrino carries away part 
of the energy, the proper decay length cannot be accu- 
rately measured. Instead, a visible proper decay length 
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(VPDL) was used in this analysis. It is defined as 



|> 30000 
O 25000 
§ 20000 
° 15000 
S 10000 
^ 5000 




3.135 0.14 0.145 0.15 0.155 0.16 0.165 0.17 

M^-WWGeV/c 2 ) 



M B (Z^.P^°)/|P£ D °| 2 - 



(3) 



Here L xy is a vector in the plane perpendicular to the 
beam direction from the primary to the B meson decay 

vertex. The transverse momentum Pt^ D is defined as 
the vector sum of the transverse momenta of the muon 
and D°. Mb is the mass of B meson. 



DESCRIPTION AND COMBINATION OF 
FLAVOR TAGGERS 



FIG. 2: The M(Kit-k) - M(Kir) invariant mass distribution 
for selected /j,D* candidates. The curve shows the result of 
the fit described in Sec lVIII 
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FIG. 3: The Kir invariant mass distribution for selected fiD° 
candidates. The curve shows the result of the fit described in 
SecEIIl 



Many different properties can be used to identify the 
initial flavor b or b of a heavy quark fragmenting into 
a reconstructed B meson. Some of them have strong 
separation power, while others are weaker. In all cases, 
their combination into a single tagging variable gives a 
significantly better result than that in the case of their 
separate use. We build such a combination with the like- 
lihood ratio method described below. We first describe 
the combination algorithm and then discuss the discrim- 
inating variables used. 



A. Combination of variables 

We construct a set of discriminating variables x\, ..., X n 
for a given event. A discriminating variable, by defini- 
tion, should have different distributions for b and b fla- 
vors. For the initial b quark, the probability density func- 
tion (p.d.f.) for a given variable Xi is denoted as flixi), 
while for the initial b quark it is denoted as f^{xi). The 
combined tagging variable r is defined as: 



n 

i=l 



(4) 
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FIG. 4: The Ktt invariant mass for selected fiD* candidates. 
The curve shows the result of the fit described in Sec lVIII 



A given variable Xi may not be defined for some events. 
For example, there are events that do not contain an 
identified muon on the opposite side. In this case, the 
corresponding variable is set to 1. An initial b flavor is 
more probable if r < 1, and a b flavor is more probable 
if r > 1. By construction, an event with r < 1 is tagged 
as a 6 quark and an event with r > 1 is tagged as a b 
quark. For an oscillation analysis, it is more convenient 
to define the tagging variable as: 



1 



1 



(5) 



By construction, the variable d which ranges between — 1 
and 1. An event with d > is tagged as a & quark and 
with d < as a b quark, with higher \d\ values corre- 
sponding to higher tagging purities. For uncorrelated 
variables x\, ...,x n , and perfect modeling of the p.d.f., d 
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gives the best possible tagging performance, and its abso- 
lute value provides a measure of the dilution of the flavor 
tagging defined in Eqn. 

Very often, the analyzed events are divided into sam- 
ples with significantly different discriminating variables 
and tagging performances. This division would imply 
making a separate analysis for each sample and combin- 
ing the results at a later stage. In contrast to this ap- 
proach, the tagging variable d defined by Eqs. 01 and 
provides a "calibration" for all events, regardless of their 
intrinsic differences. Since the absolute value of d gives a 
measure of the dilution of the flavor tagging, events from 
different categories but with a similar absolute value of d 
can be treated in the same way. Thus, another important 
advantage of this method of flavor tagging is the possibil- 
ity of building a single variable having the same meaning 
for different kinds of events. It allows us to classify all 
events according to their tagging characteristics and use 
them simultaneously in the analysis. 

All of the discriminating variables used in this anal- 
ysis are constructed using the properties of the b quark 
opposite to the reconstructed B meson ("opposite-side 
tagging"). Since an important property of the opposite- 
side tagging is the independence of its performance of 
the type of reconstructed B meson, it can be calibrated 
in data by applying tagging to the events with B° and 
B + decays. The measured performance can then be used 
to study B® meson oscillations, as an example. 

The probability density functions for each discrimi- 
nating variable discussed below were constructed using 
events from the D° sample with < VPDL < 500 /im. 
In this sample, the decay B + — > /.i + vD° dominates, see 
Sec. rvTTH The B Q d -► p+vD*+ events give a 16% con- 
tribution to the sample and, due to the cut on VPDL, 
contains mainly non-oscillated B® decays, as determined 
by Monte Carlo (the standard pythia generation, fol- 
lowed by decay of B mesons with EvtGen, passed through 
Geant and then reconstruction). 

The initial flavor of a b quark is therefore determined 
by the charge of the muon. Estimates based on Monte 
Carlo simulation indicate that the purity of the initial fla- 
vor determination in the selected sample is 0.98 ± 0.01, 
where the uncertainty is due to the uncertainties in mea- 
sured branching fractions of B meson decays. 

For each discriminating variable, the signal band 
containing all events with 1.80 < M(Kn) < 1.92 
GeV/c 2 and the background band containing all events 
with 1.94 < M(Ktt) < 2.2 GeV/c 2 were defined. The 
p.d.f's were constructed as the difference in the distribu- 
tions. The latter distributions were normalized by multi- 
plying them by 0.74 so that the number of events in the 
background band corresponds to the estimated number 
of background events in the signal band. 



B. Flavor tagger discriminants 

We now describe the variables used. An additional 
muon was searched for in each analyzed event. This muon 
was required to have at least one hit in the muon cham- 
bers and to have cos 4>(p^,pb) < 0.8, where ps is the 
three-momentum of the reconstructed B meson, and <f> is 
the angle between the vectors p M and ps- If more than 
one muon was found, the muon with the highest number 
of hits in the muon chambers was used. If more than 
one muon with the same number of hits in the muon 
chambers was found, the muon with the highest trans- 
verse momentum px was used. For this muon, a muon 
jet charge Qj was constructed as 



Qj 



where q l is the charge and p l T is the transverse momentum 
of the i'th particle, and the sum is taken over all charged 
particles, including the muon, satisfying the condition 
AR = ^{A(f>) 2 + (An) 2 < 0.5, where A<j) and An are 
computed with respect to the muon direction. Daughters 
of the reconstructed B meson were explicitly excluded 
from this sum. In addition, any charged particle with 
cos cf>(p, pb) > 0.8 was excluded. The distribution of 
the muon jet charge variable is shown in Figs. Ela) and 
(b). In these plots, q roc gives the charge of the b quark 
in the reconstructed B — > p + vD° decay, in this case 
given by the muon charge. We build separate p.d.f.'s for 
muons with hits in all three layers of the muon detector, 
Fig. Ufa), and for muons with fewer than three hits, Fig. 

In addition to the muon tag, reconstructed electrons 
with cos 4>(p e ,pB) < 0.8 were also used for flavor tag- 
ging. The electron is reconstructed by extrapolating a 
track to the calorimeter and adding up the energy de- 
posited in a narrow tube or "road" around the track. 
Calorimeter cells are collected around the track extrap- 
olated positions in each layer and the total transverse 
energy of the cluster is defined by the sum of the en- 
ergies in each layer. The electrons are required to be 
in the central region (\n\ < 1.1), with px > 2 GeV/c. 
They are required to have at least one hit each in the 
CFT and SMT. They are required to have energy de- 
posits in the EM calorimeter consistent with an electron, 
0.55(0.5) < E/p < 1.0(1.1), and low energy deposit in 
the hadron calorimeter, EMF > 0.8(0.7). The cuts are 
looser for electrons with px > 3.5 GeV/c and are given 
in brackets. EMF and E/p are calculated as below: 



EMF = 



E/p 



layer 



layer number i— 1,2,3 



E T {i) 



Salllayers-^rCO 
^^layer number i— 1,2,3 -^t(^) 



px(track) 



(6) 
(7) 



where Et(i) is the transverse energy within the road 
in the i'th layer. We also require a minimum single 
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layer cluster energy of a cluster in the central preshower, 
CPS| LC > 4.0(2.0) MeV/c. The cuts were optimized by 
studying electrons from conversion decays (7 — > e + e~) 
and fakes from Kg — > ir + Tr~ decays to obtain a 90% 
purity for electrons. For these electrons, an electron jet 
charge (Qj) was constructed in the same way as the muon 
jet charge, Q^j. The distribution of the electron jet charge 
variable is shown in Fig.lSTc). 
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FIG. 5: (a) Distribution of the jet charge Qj for muons with 
hits in all three layers of the muon detector, (b) Distribution 
of the jet charge Q j, for muons with fewer than three hits, (c) 
Distribution of the jet charge for electrons Qj. Here q(b rcc ) 
is the charge of the muon from the reconstruction side. 



secondary vertex identification algorithm can be found in 
Rcf. 0. 

The three-momentum of the secondary vertex psv is 
defined as the vector sum of the momenta of all tracks 
included in the secondary vertex. A secondary vertex 
with cos 4*(psVtPb) < 0.8 was used for flavor tagging. 
A secondary vertex charge Qsv is defined as the third 
discriminating variable 



Qsv 



where the sum is taken over all tracks included in the sec- 
ondary vertex. Daughters of the reconstructed B meson 
were explicitly excluded from this sum. In addition, any 
charged particle with cos0(p,ps) > 0.8 was excluded. 
Here p L is the longitudinal momentum of track i with re- 
spect to the direction of the secondary vertex momentum 
Pv- A value of k = 0.6 was used, taken from previous 
studies at LEP Q ■ We verified that this value of k results 
in the optimal performance of the Qsv variable. Figures 
HJa) and H^b) show the distribution of this variable for 
the events with and without an identified muon flavor 
tag. 

Finally, the event charge Qev was constructed as 



Qev 



Si i 1 pt 



The sum is taken over all charged tracks with 0.5 < px < 
50 GeV/c and having cos 4>(p,pb) < 0.8. Daughters of 
the reconstructed B meson were explicitly excluded from 
this sum. The distribution of this variable is shown in 
Fig. Etc). 

For each event with an identified muon, the muon jet 
charge Qj and the secondary vertex charge Qsv were 
used to construct a muon tagger. For each event without 
a muon but with an identified electron, the electron jet 
charge Qj and the secondary vertex charge Qsv were 
used to construct an electron tagger. Finally, for events 
without a muon or an electron but with a reconstructed 
secondary vertex, the secondary vertex charge Qsv an d 
the event jet charge Qev were used to construct a sec- 
ondary vertex tagger. The resulting distribution of the 
tagging variable d for the combination of all three tag- 
gers, called the combined tagger, is shown in Fig. 
The performances of these taggers are discussed in the 
following sections. 



An additional secondary vertex corresponding to the 
decay of a B hadron was searched for, using all charged 
tracks in the event excluding those from the recon- 
structed B hadron. The secondary vertex was also re- 
quired to contain at least two tracks with an axial im- 
pact parameter significance greater than 3. The distance 
l xy from the primary to the secondary vertex must also 
satisfy the condition: l xy > ia(l xy ). The details of the 



VI. MULTIDIMENSIONAL TAGGER 

In addition to the flavor tagger described in Sec. [V] 
an alternative algorithm was also developed and used to 
measure B° mixing. This tagger is multidimensional, 
i.e., the likelihood functions it is based on depend on 
more than a single variable. In addition, the p.d.f.'s were 
determined from simulated events, while the primary fla- 
vor tagger described in Sec[V]uses data to construct the 
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FIG. 6: (a) Distribution of the secondary vertex charge for 
events with an opposite-side muon. (b) Distribution of the 
secondary vertex charge for events without an opposite-side 
muon. (c) Distribution of the event jet charge. q{b IOC ) is the 
charge of the b quark from the reconstruction side. 



p.d.f.'s. The multidimensional tagger therefore provides 
a cross-check of the primary algorithm. 

If, as before, we have a set of discriminants x%, ...,x ni 
the likelihood that the meson has flavor b at the time 
of creation can be written as C(b] xi, x n ). A similar 
expression C(b; X\, x n ) holds for the likelihood for b. 
These likelihoods relate to the variable d as 



d 



C(b)+C(b)' 



(8) 



This definition is similar to Eq. JJjJ. 

The likelihoods arc obtained from the simulated sam- 
ples of B ± -> J/tfjK* 1 with J/iji -> ^V~- This final 
state does not oscillate and is therefore flavor-pure. The 
B~ — > J/^K~ sample was used to obtain £(&), while 
C(b) was determined from B + — > J/%pK + sample. In 
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FIG. 7: Normalized distributions of the combined tagging 
variable. q(b lec ) is the charge of the 6 quark from the recon- 
struction side. 



practice, the likelihoods were stored as multidimensional 
histograms (with one dimension per discriminating vari- 
able) with the bin content normalized to the total num- 
ber of events in the sample. For a given event, the tagger 
output d was obtained by substituting the appropriate 
normalized bin contents into Eq. JSJ. 

In addition to the discriminating variables introduced 
in Sec0 other variables were used for the multidimen- 
sional tagger. For each identified opposite-side muon, 
the transverse momentum px relative to the beam axis 
and transverse momentum p^ 1 relative to the nearest jet 
were computed. (The muon was included in the jet clus- 
tering.) Another variable defined for the muon is its im- 
pact parameter significance S^, where 5* M is the trans- 
verse impact parameter significance eT/c(eT), where ex 
is defined in Sec. IIHI For each reconstructed opposite- 
side secondary vertex, the secondary vertex transverse 
momentum p T v was computed by taking the magnitude 
of the transverse projection of the vector sum of all tracks 
in that vertex. In principle, all discriminating variables 
can be combined into a single multidimensional likeli- 
hood. However, since a binned likelihood was used, in 
order to achieve a reasonable resolution in any given dis- 
criminant, the binning must be fine enough to resolve its 
useful features. In practice, because of limited simulation 
statistics, this means that discriminating variables must 
be chosen wisely when making a combination. 

All events were divided into three categories based on 
their opposite-side content. The following variables for 
different categories were selected. 



1. Events with muon and secondary vertex: 
Tag(fi+SV)={Q^pf;Q sv }. 

2. Events with muon and without secondary vertex. 
Tag( Ai -SV)={Q^;^ 1 ;p T ;^ I }. 

3. Events with secondary vertex without a muon: 
Tag(SV-n)={Q E v,QsV,P% v }- 

Distributions in the tagging variable d for the above 
three taggers are shown in Fig. |SJ They were made by 
applying the taggers to the simulated B ± — > J/tpK^ 
samples from which they were created. 

The final multidimensional tagger used the following 
logic to decide which of its sub-taggers to use. For 
events containing a muon and a secondary vertex, the 
Tag(/x+SV) was used. If the opposite side contained 
a muon and no secondary vertex, the Tag(/i— SV) was 
used. If the opposite side contained an electron, the elec- 
tron tagger described in Sec. IVlwas used. Note that this 
tagger is not multidimensional and is not derived from 
simulation. If the opposite side contained a secondary 
vertex, the Tag(SV— /i) was used. 



VII. ASYMMETRY FIT PROCEDURE 

The performance of the flavor tagging and measure- 
ments of the B° mixing frequency Am^ were obtained 
from a study of the dependence of the flavor asymmetry 
on the B-meson decay length. 

The flavor asymmetry A is defined as: 

_ A nos - N osc 

Here N nos is the number of non-oscillated B decays and 
iV osc is the number of oscillated B decays. An event 
B — > fj, + vD° X with x d < was tagged as non- 
oscillated, and an event with x d > was tagged as 
oscillated. The flavor tagging variable d is defined in Eq. 
© or ©. 

All events in the D° and D* samples were divided into 
seven groups according to the measured VPDL (L) de- 
fined in Eq. ©. The numbers of oscillated JV° SC and 
non-oscillated N^° s signal events in each group i were 
determined from the number of the D° signal events 
given by a fit to the Kir invariant mass distribution for 
both samples. The seven VPDL bins (in cm) defined 
were: 

-0.025 < L < 0.0, 0.0 < L < 0.025, 0.025 < L < 0.050, 
0.050 < L < 0.075, 0.075 < L < 0.1, 0.1 < L < 0.125 
and 0.125 < L < 0.2. 



A. Mass Fit 

In this section we describe the mass fitting procedure. 
The fitting function was chosen to give the best x 2 °f 
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FIG. 8: Normalized distributions of the combined tagging 
variable for three multidimensional taggers for the simulated 
samples B — > J/jpK ± . Here q(b TCC ) is the charge of the b 
quark from the reconstructed side, (a) Distribution of d for 
Tag(^+SV). (b) Distribution of d for Tag(^-SV). (c) Distri- 
bution of d for Tag(SV-^i). 



the fit to the A'tt mass spectrum of the entire sample of 
B — > fi + D° X events shown in Figs. |3]and0] The signal 
peak corresponding to the decay D° — » K~ix + can be 
seen at 1.857 GeV/c 2 . The background to the right of the 
signal region is adequately described by an exponential 
function: 

f**{x) = a x e -S , (10) 

where x is the Kir mass. 

The peak in the background to the left of the signal is 
due to events in which D mesons decay to KttX where 
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X is not reconstructed. It was modeled with a bifurcated 
Gaussian function: 



/ 2 bkg M = A e 



A e 



(ic— /to) 



(2(T 



,) 



for x- no > (11) 
for x — fig < 0. 



Here /io is the mean of the Gaussian, and <jl and cxr are 
the two widths of the bifurcated Gaussian function. 

The signal has been modeled by the sum of two Gaus- 
sians: 



f is (2) 



iV slg / n 

\72vr \ cri ' 



0"2 



,(12) 



where N slg is the number of signal events, fj-i an d A*2 are 
the means of the Gaussians, <j\ and <ii are the widths 
of the Gaussians, and r\ is the fractional contribution of 
the hrst Gaussian. 

The complete fitting function, which has twelve free 
parameters, is: 



(13) 



The low statistics in some VPDL bins, which have as 
few as ten events after flavor tagging, do not permit a 
free fit to this function, Consequently some parameters 
had to be constrained or fixed. In order to do this, it was 
necessary to show that the constraints on the parameters 
are valid for all of the VPDL bins. Unconstrained fits 
were performed to several high statistic samples, and the 
set of all events was used as a reference fit. Events were 
divided into VPDL bins and fit to investigate the VPDL 
dependance of the fit results. In addition, three samples 
were made to test whether the presence of a flavor tag 
changes the mass spectrum: all tagged events over the 
entire VPDL range, all events in the short VPDL range 
[0,0.05] tagged as opposite-sign events, and all events in 
VPDL range [0,0.05] tagged as same-sign events. 

This study showed that the width, position, and the 
ratio of the signal Gaussians, as well as the position and 
widths of the bifurcated Gaussian describing the back- 
ground can be fixed to the values obtained from the fit 
to the total D° or D* mass distribution. This left four 
free parameters: the numbers of events in the signal peak, 
background peak, and exponential background, and the 
slope constant of the exponential background. Examples 
of the fits to the Kir mass distribution in different VPDL 
bins are shown in Fig. [5] 

The number of D* candidates was estimated using the 
distribution of (Mk+^-k- —Mr+„-), shown in Fig. [5] In 
this case, the signal was modeled with two Gaussians as 
described by Eci. (|12fl . and the background by the product 
of a linear and exponential function 

/ bkg (:r) = a[l + c(x-x )]e^. (14) 

where x is the mass difference {M K +^- v - — Mj^+^~) in 
this equation. 
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FIG. 9: The fit to M(K-k) mass for non-oscillating (left) and 
oscillating (right) B — > fi + i>D*~ events tagged by the muon 
tagger with \d\ > 0.3 in VPDL bins 0.0 < L < 0.025 cm 
(2a,2b), 0.025 < L < 0.050 cm (3a,3b), 0.075 < L < 0.100 cm 
(5a,5b), and, 0.100 < L < 0.125 cm (6a,6b). 



B. Expected flavor asymmetry 



For a given type of B q meson (q — u,d,s), the dis- 
tribution of the visible proper decay length L is given 
by 
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<° S {L,K) 

n° d sc (L,K) 

nT S (L,K) 
n° s sc (L,K) 



K 



1 

2 ' ct{B+)' 
1 K 



^J(l+V u ), (15) 



2 ct(B+) 

i is: 



^>){X-V U ), (16) 



2 ct( j B°) 

I£ L , * 

1 +X>d cos(Am d ) 

c 



1 



A" 



2 ct(B°) 
1 — 2?d cos(Arn,i 



A A 



1 A 

1 A" 



2 CT (B0) 



(17) 

(18) 
(19) 
(20) 



Here r is the lifetime of the B meson, Am^ is the mix- 
ing frequency of B° mesons, the factor A = PjP / 

reflects the difference between the measured (P£ D ) and 
true (Bp) momenta of the B meson. The B + meson 
does not oscillate, and it is assumed in these studies that 
the B® meson oscillates with infinite frequency. The fla- 
vor tagging dilution is given by T>. In general, it can be 
different for B° and B + . In our study we verified the 
assumption that T>d = T> u for our opposite-side flavor 
tagging. 

The transition from the true to the experimentally 
measured visible proper decay length L is achieved by 
integration over the A-factor distribution and convolu- 
tion with the resolution function: 



nos/osc^M) 



= J dL Rj(L - L M ) Ej (L) 9(L) 
x J dK Dj(K) n n q °) /osc (L,K) (21) 



Here Rj(L — L ) is the detector resolution in the 
VPDL, and Ej(L) is the reconstruction efficiency for a 
given channel j of B q meson decay. The step function 
9(L) forces L to be positive in the integration. L M can be 
negative due to resolution effects. The function Dj(K) 
is a normalized distribution of the A"-factor in a given 
channel j, obtained from simulated events. 

In addition to the main decay channel B — > /i + ^£> X, 
the process cc — > n + vD Q X contributes to the selected 
final state. A dedicated analysis was developed to study 
this process, both in data and in simulation. It shows 
that the pseudo decay-length, constructed from the cross- 
ing of the \i and D trajectories, is distributed around 
zero with a w 150 /im. The distribution N c£ (L M ) of the 
VPDL for this process was taken from simulation. It was 



assumed that the production ratio (c — > D*)/(c — > D°) is 
the same as in semileptonic B decays and that the flavor 
tagging for the cc events gives the same rate of oscillated 
and non-oscillated events. The fraction f cS of cc events 
was obtained from the fit. 

Taking into account all of the above mentioned contri- 
butions, the expected number of (non-) oscillated events 
in the i-th bin of VPDL is 



N 



e.nos/osc 



dL M (1 - f a 



yq—u,d.s j 

- M 



N 



+ I dL M f c5 N c5 (L M ). 



nos/osc 



(22) 



Here the integration J. dL M is taken over a given interval 
i, the sum is taken over all decay channels B q — > 

fj, + vD° X contributing to the selected sample, and Br,- is 
the branching fraction of channel j. 

Finally, the expected value of asymmetry, 
vl-(Am, f c5 , Dd, T> u ), for the interval i of the mea- 
sured VPDL is given by 



At(AmJ cS ,V d ,V u ) 



NT 



Nf 



NT 



NT 



(23) 



The expected asymmetry can be computed both for 
the D* and the D° samples. The only difference between 
them is due to the different relative contributions of var- 
ious decay channels of B mesons. 

For the computation of Af , the B meson lifetimes and 
the branching fractions Br, were taken from the Particle 
Data Group (PDG) |3J. They are discussed in the follow- 
ing section. The functions Dj(K), Rj(L), and Sj(L) were 
obtained from MC simulation. Variations of these inputs 
within their uncertainties are included in the systematic 
uncertainties. 



C. Sample Composition 

There is a cross-contamination between the b — > B° — > 
H+vD°X, b -► B° -► /i+i/j) !; and b -> B+ 
fj, + vD° X samples. To determine the composition of the 
selected samples, we studied all possible decay chains for 
-B , B®, and B + with their corresponding branching frac- 
tions, from which we estimated the sample composition 
in the D* and D° samples. 

The following decay channels of B mesons were con- 
sidered for the D* sample: 



B a 
B° 
B+ 
B° 



fi+ivD** - 



X: 
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and for the D° sample: 



B+ - 






B+ - 


■* v+vD* , 




B+ - 


■* fi+vD**° - 


> /x+zaD X, 


B+ - 


-» fi + uD**°- 


► fj, + vD*°X, 


B° - 


-> n + vD**- - 


-> /x+zaD° X, 


B° - 


-> (i+vD**- - 


-» n+^D*°X, 


B° s - 


-» fi+i>D X, 




B° s - 


-> li + vD*°X. 





Here, and in the following, the symbol "£)**" denotes 
both narrow and wide D** resonances, as well as non- 
resonant -D-7T and _D* 7r production. 

The most recent PDG values were used to determine 
the branching fractions of decays contributing to the I? 
and -D* samples: 

Br(B+ -> /i+i/.D Q ) = (2.15 ± 0.22) x 10" 2 , 

Br(_B° — > n + vD~) = (2.14 ± 0.20) x 10" 2 , 

Br(f? + — > i± + vD*°) = (6.5 ±0.5) x 10~ 2 , 

Br(B° -> n + vD*-) = (5.44 ± 0.23) x 10" 2 . 

Br(£> + — > n + vD**°) was estimated using the following 
inputs: 

Br(B -> /x+iaX") = (10.73 ± 0.28) x 10" 2 0, 
Br(B°^/i+^) = t{B°)/t{B + ) 

■ Br(B+ -> jU+2/X), 
Br(B+ -> n + vD**°) = Br(S+ -> ,u + z/X) 

- Br(B+ -> n + vD°) 

- Br(S+ -> n+uD*°), (24) 

where t(B°) is the -B lifetime, and r(_B + ) is the -B + 
lifetime. The following value was obtained: 

Br(B+ -► fi+vD** ) = (2.70 ± 0.47) x 10~ 2 . (25) 

Br(i?° — > n + vD**~) is obtained as follows: 

Br(B° ij+vD**-) = Br(B° -> /x+i/X) 

- Br(B° -> n + vD~) 

- Br(B 
r(B° 



Br(£r -> fj, + i/D**~ 

Br(S+ -> f i+vD**° 
the following inputs: 



7!/ ,. Br(B + — * /x + z/.D ) 
vD*~ X) was estimated from 



Br(6 -> l + vD*~ ti + X) = (4.73 ±0.8 ±0.6) x 10" 3 [TJ, 
Br(6 -> l + vD*- ti + X) = (4.80 ±0.9 ±0.5) x 10" 3 H3, 
Br(5 -> 1+vD*-it-X) = (0.6 ±0.7 ±0.2) x 10" 3 [13 

and assuming Br(6 -> B+) = (0.397 ± 0.010) |]. The 
usual practice in estimating this decay rate is to neglect 



the contributions of the decays D** — > D*irir. However, 
the above data allows us to take these decays into ac- 
count. 

Neglecting the decays D** — > D*iririr, the available 
measurements can be expressed as: 



Br(S -> 1 + vD*-it + X) 



Br(B+ 
Br(5° 



► l + vD*-n+X a ) 

l + l>D*-Tt + TT-). 



Br{B -> I+vD^tt^X) = Br(S° -> 

From these relations and using the above measurements, 
we obtain 

Br(B+ -► ii+vD**° -> l+vD*-X) = 

(1.06 ± 0.24) x 10" 2 . (26) 

All_ other factors for the Br (5 -> fj+uD** -> 
fj, + vD* X) were obtained assuming the following rela- 
tions, 



Br(S° -> /i+i/L>* 



j.i + vD* 7r) 



Br(B+ ->■ n+vD**° -> ^i+i/D* 7r) 
Br(£ -► ^+i/L>** -► ^+i/L>* 7T+) 



= r{S°)/r(B+) 1 
= 2. 



Br(B -> n+vD** -> tt°) 

fi + vDX) was estimated from the 



Br(B -> n+vD** 
following inputs: 



Br(B -> [i+vD** -> ii+vDtt+) 
Br(B {l+vD** ii+vD-k ) 



= 2, 



Br(B -> [i+vD**) = Br(B ^i+vD** ^i+vDX) 
+ Br(B ii + vD** ii + vD* X). 

To estimate branching fractions for decays, Br(£?° — > 
Li+vD-X) = (7.9 ± 2.4) x 10" 2 was taken from Ref. 
and the following assumptions were used: 



Br(B° 
Br(B s ° 



Br(B° n+vDT' 



H+vD*°X) ' 



where r(B°) is the B° meson lifetime. In addition, it was 
assumed that 

Br^g^+^r-^+^P**) _ n « , 97A 
Br(flO _ n+yD?) ~ [ ' 

There is no experimental measurement of this ratio yet 
and to estimate the the corresponding systematic uncer- 
tainty, this ratio was varied between and 1. 

In addition to these branching fractions, various de- 
cay chains are affected differently by the B meson selec- 
tion cuts, and the corresponding reconstruction efficien- 
cies were determined from simulation to correct for this 
effect. Taking into account these efficiencies, the compo- 
sition of the D* sample was estimated to be (0.89±0.03) 
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B°, (0.10 ± 0.03) B+, and (0.01 ± 0.01) B° s . The D° 
sample contains (0.83 ± 0.03) B+ , (0.16 ±0.04) B° , and 
(0.01 ±0.01) B°. 

Since the D* sample was selected by the cut on the 
mass difference AM = M(D°ir) - M(D°), there is a 
small additional contribution of B — > n + uD° events to 
the D* sample, when D° is randomly combined with a 
pion from the combinatorial background. The fraction 
of this contribution was estimated using fj, + D° ir + events. 
These events were selected applying all the criteria for the 
D* sample, described in Sec. IIIII except that the wrong 
charge correlation of muon and pion was required, i.e., 
the muon and the pion were required to be having the 
same charge. The number of D° events was determined 
using the same fitting procedure as for the D* sample, 
and the additional fraction of B — > ^i + vD° events in 
the D* sample was estimated to be (4.00 ± 0.85) x 10~ 2 . 
This fraction was included in the fitting procedure and 
the uncertainty in this value was taken into account in 
the overall systematics. 



VIII. RESULTS 

For each sample of tagged events, the observed and 
expected asymmetries were determined using Eqs. (|5J) 
and (|23[1 in all VPDL bins, and the values of Am d , fee, 
T> u , and T> d were obtained from a simultaneous x 2 fit: 

X 2 (Am d , fcc,'D d ,T> u ) = 

X 2 D * (Am d , f cS ,T> d ,T> u ) +Xl>o (Amd, fee, T> d , V u )\ (28) 
X%* {Am d J C c,V d ,'D u ) = 

\ A i,D* - A\ D , {Am d J cS ,V d ,V u )} 2 



E 



X 2 D o (Am d , f cS ,V d ,V u ) 



E 



[A,n° -A e iD0 (Am* f c - c ,V d ,V u )f 



(29) 



(30) 



Here is the sum over all VPDL bins. Examples of the 
X 2 fit to the flavor asymmetry minimization given in Eq. 
(|2~5jl . is shown in Fig. HTH 

The performance of the flavor tagging method was 
studied separately for the muon, electron, and secondary 
vertex taggers using events with \d\ > 0.3. Results are 
given in Tables HI IIIII All uncertainties in these tables are 
statistical only and do not include systematic uncertain- 
ties. The performances of the combined tagger defined 
in Sec lVBI for events with \d\ > 0.3 and the alternative 
multidimensional tagger defined in Sec I VII for events with 
\d\ > 0.37 are also shown. The cut on |d| is somewhat 
different for the multidimensional tagger as the calibra- 
tion is different and we compare the dilutions for the same 
tag efficiency for the two taggers. The tagging efficiencies 
shown in TablesHlandlTTIwere computed using events with 
VPDL=[0. 025, 0.250]. This selection reduces the contri- 
bution from cc — > ^, + vD°X events, since they have a 



VPDL distribution with zero mean and a w 150 /jm as 
described in Sec lVllI 

Individual taggers give compatible values of Am d and 
/ c g, as can be seen in Table ITTT1 For the combined tagger 
with |d| > 0.3, the following results were obtained: 



eV\ = (2.19 ±0.22)%, 
Am d = 0.513 ± 0.023 ps~ 
fee = (3.3 ±1.3)%. 



(31) 



The multidimensional tagger which used simulation for 
the description of p.d.f.'s as described in Sec I VII gives 
consistent results both for the Am d and the fraction 
f cS , which is used as a cross-check of the main tagging 
algorithm. 

One of the goals of this analysis was to verify the as- 
sumption of independence of the opposite-side flavor tag- 
ging on the type of the reconstructed B meson. It can 
be seen from Tables Q] and [n] that the measured flavor 
tagging performance for B° events is slightly better than 
for B + events, both for individual and combined taggers. 
This difference can be explained by the better selection 
of fi + i/D*~ events due to an additional requirement of 
the charge correlation between the muon and pion from 
D*~ — ► D° 7r~ decay. The D° sample can contain 
events with a wrongly selected muon. Since the charge 
of the muon determines the flavor asymmetry, such a 
background can reduce the measured B + dilution. The 
charge correlation between the muon and the pion sup- 
presses this background and results in a better measure- 
ment of the tagging performance. 

To verify this hypothesis, a special sample of events 
satisfying all conditions for the D* sample, except the re- 
quirement of the charge correlation between the muon 
and the pion, was selected. The dilution V' d for this 
sample is shown in Table [H] It can be seen that T>' d 
is systematically lower than T> d for all samples and all 
taggers. T>' d is the right quantity to be compared with 
T> u and Table ITT1 shows that they are statistically com- 
patible. This result therefore confirms the expectation of 
the same performance of the opposite-side flavor tagging 
for B + and B° events. It also shows that contribution of 
background in the D° sample reduces the measured di- 
lution for B + events. Thus, the dilution measured in the 
D* sample can be used for the B® mixing measurement, 
where a similar charge correlation between the muon and 
D s is required. 

By construction, the dilution for each event should 
strongly depend on the magnitude of the tagging variable 
d. This property becomes important in the B® mixing 
measurement, since in this case the dilution of each event 
can be estimated using the value of d and can be included 
in a likelihood function, improving the sensitivity of the 
measurement. To test the dependence of the dilution on 
d, all tagged events were divided into subsamples with 
0.1 < |d| < 0.2, 0.2 < \d\ < 0.35, 0.35 < \d\ < 0.45, 
0.45 < \d\ < 0.6, and |d| > 0.6. The overall tagging effi- 
ciency for this sample is (19.95 ± 0.21)%. The dilutions 
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FIG. 10: Asymmetries obtained in the D* and D° sample with the combined tagger in |d| bins. Circles are data, and the result 
of the fit is superimposed. 
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TABLE I: Tagging performance for the D* sample for different taggers and subsamples. Uncertainties are statistical only. 



Tagger 


e{%) 


v d 


eV 2 d (%) 


Muon (\d\ > 0.3) 


6.61 ±0.12 


0.473 ± 0.027 


1.48 ±0.17 


Electron (\d\ > 0.3) 


1.83 ±0.07 


0.341 ± 0.058 


0.21 ±0.07 


SVCharge (\d\ > 0.3) 


2.77 ±0.08 


0.424 ± 0.048 


0.50 ±0.11 


Combined (\d\ > 0.3) 


11.14 ±0.15 


0.443 ± 0.022 


2.19 ±0.22 


TV T li ■ 1 ■ / 1 J 1 ^ n (1I7\ 

Multidim. (\d\ > 0.37) 


10.98 ± 0.15 


0.395 ± 0.022 


1.71 ± 0.19 


Combined (0.10< \d\ <0.20) 


4.63 ±0.10 


0.084 ± 0.031 


0.03 ±0.02 


Combined (0.20< \d\ <0.35) 


5.94 ±0.12 


0.236 ± 0.027 


0.33 ±0.08 


Combined (0.35< \d\ <0.45) 


3.89 ± 0.09 


0.385 ± 0.034 


0.58 ±0.10 


Combined (0.45< \d\ <0.60) 


4.36 ±0.10 


0.512 ±0.032 


1.14 ±0.14 


Combined (0.60< \d\ <1.00) 


1.13 ±0.05 


0.597 ± 0.058 


0.40 ± 0.08 



TABLE If: Tagging performance for the D° sample for different taggers and subsamples. For comparison, the dilution T>' d 
measured in the D* sample with addition of wrong sign /j, + uD° tv + events is also shown. Uncertainties are statistical only. 



Tagger 


s{%) 


T> u 


eVl{%) 


v' d 


Muon (\d\ > 0.3) 


7.10 ±0.09 


0.444 ± 0.015 


1.400 ±0.096 


0.463 ± 0.028 


Electron (\d\ > 0.3) 


1.88 ±0.05 


0.445 ± 0.032 


0.372 ± 0.054 


0.324 ± 0.060 


SVCharge (|d| > 0.3) 


2.81 ±0.06 


0.338 ± 0.026 


0.320 ± 0.050 


0.421 ± 0.049 


Combined > 0.3) 


11.74 ±0.11 


0.419 ±0.012 


2.058 ±0.121 


0.434 ± 0.023 


Multidim. (\d\ > 0.37) 


11.67±0.11 


0.363 ±0.012 


1.540 ±0.106 


0.384 ± 0.023 


Combined (0.10< \d\ <0.20) 


4.59 ± 0.08 


0.104 ±0.017 


0.050 ±0.016 


0.079 ± 0.029 


Combined (0.20< \d\ <0.35) 


6.10 ±0.09 


0.234 ±0.014 


0.335 ± 0.042 


0.212 ±0.024 


Combined (0.35< \d\ <0.45) 


3.98 ± 0.07 


0.361 ±0.018 


0.519 ±0.052 


0.364 ± 0.032 


Combined (0.45< \d\ <0.60) 


4.77 ±0.07 


0.504 ±0.016 


1.211 ±0.077 


0.489 ± 0.030 


Combined (0.60< \d\ <1.00) 


1.17 ±0.04 


0.498 ± 0.031 


0.290 ±0.038 


0.572 ± 0.056 



TABLE III: Measured value of Am d and / C c for different 
taggers and subsamples. 



Tagger 


Am d 


(ps 


- 1 ) 




fee 






Muon 






0.502 


±0 


028 





013 


± 





010 


Electron 






0.481 


±0 


067 





058 


± 





045 


SVCharge 






0.553 


±0 


053 





096 


± 





050 


Multidim. 






0.502 


±0 


026 





031 


± 





014 


Combined 


> 0.3) 




0.513 


±0 


023 





033 


± 





013 


Combined 


(0.10 < \d 


< 0.20) 


0.506 


±0 


209 





495 


± 





505 


Combined 


(0.20 < \d 


< 0.35) 


0.523 


±0 


064 





021 


± 





025 


Combined 


(0.35 < \d 


< 0.45) 


0.531 


±0 


042 





063 


± 





038 


Combined 


(0.45 < \d 


< 0.60) 


0.510 


±0 


032 





010 


± 





010 


Combined 


(0.60 < \d 


< 1.00) 


0.456 


±0 


049 





032 


± 





026 



obtained are shown in Table [H Their strong dependence 
on the value of the tagging variable is clearly seen. This 
allows us to perform a dilution calibration and obtain the 
measured dilution T>d as a function of the predicted value 
This is used to provide an event-by-event dilution for 
the B s mixing analysis and the calibration derived in this 
analysis is used for the two sided C.L. on B s mixing, ob- 



tained by D0 . The overall tagging power, computed 
as the sum of the tagging powers in all subsamples, is: 

eV\ = (2.48 ±0.21)%. (32) 

The measured oscillation parameters Am^ for all 
considered taggers and subsamples are given in Ta- 
ble lllll They are compatible with the world average value 
Am d = 0.502 ± 0.007 ps -1 3 in each instance. 

The final mixing parameter Amj was obtained from 
the simultaneous fit of the flavor asymmetry in the vari- 
ous tagging variable subsamples defined above. The frac- 
tion f cS was constrained to be the same for all subsam- 
ples. The result is 

Am d = 0.506 ± 0.020 ps" 1 (33) 

fee = (2.2 ±0.9)%. 

The statistical precision of Amd from the simultaneous 
fit is about 10% better than that from the fit of events 
with \d\ > 0.3. This improvement is directly related to 
a better overall tagging power [Eq. Q32f) ] for the sum of 
subsamples as compared to the result [Eq. (|31|l ] for the 
sample with \d\ > 0.3. 
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IX. SYSTEMATIC UNCERTAINTIES 

The systematic uncertainties are summarized in Ta- 
bles IIVI and [V] Table IIVI shows the contributions to the 
systematic uncertainty in Am^. Table Ivl shows the cor- 
responding contributions to the systematic uncertainties 
inV(B Q ). 

These uncertainties were obtained as follows: 

• The B meson branching fractions and lifetimes used 
in the fit of the asymmetry were taken from Ref. 
Q and were varied by one standard deviation. 

• The VPDL resolution obtained in simulation was 
multiplied by factors of 0.8 and 1.2. These factors 
exceed the uncertainty in the difference of the res- 
olution between data and simulation. 

• The variation of iVfactors with the change in the 
B momentum was neglected in this analysis. To 
check the impact of this assumption on the final 
result, the computation of i\-factors, was repeated 
without the cut on pt(D° ) or by applying an ad- 
ditional cut on the pt of muon, px > 4 GeV/c. 
The change in the average values of the i^-factors 
did not exceed 2%, which was used as the estimate 
of the systematic uncertainty in their values. This 
uncertainty was propagated into the variation of 
Arrid and tagging purity by repeating the fit with 
the if-factor distributions shifted by 2%. 

• The ratio of the reconstruction efficiencies in differ- 
ent B meson decay channels depends only on the 
kinematic properties of corresponding decays and 
can therefore be reliably estimated in the simula- 
tion. The ISGW2 model [3 of semileptonic B de- 
cays was used. The uncertainty in the reconstruc- 
tion efficiency, set at 12%, was estimated by vary- 
ing the kinematic cuts on the px of the muon and 
D° in a wide range. Changing the model describ- 
ing semileptonic B decay from ISGW2 to a HQET- 
motivated model [lOj produces a smaller variation. 
The fit to the asymmetry was repeated with the 
efficiencies to reconstruct the B — > n + vD**~ and 
B — > n+vD** channels modified by ±12%, and the 
difference was taken as the systematic uncertainty 
from this source. 

• The additional fraction of D° events contributing 
to the D* sample was estimated at (4.00 ± 0.85)% 
(see Sec. IVII CI) . This variation was used to esti- 
mate the systematic uncertainty from this source. 
As a cross-check, the number of D* events was 
determined from the fit of the mass difference 



M(D° ir)-M(D° ) and the fit of the flavor asymme- 
try was repeated. The measured value of Am^ = 
0.507 ± 0.020 ps" 1 is consistent with Eq. (33J. 

• We also investigated the systematic uncertainty in 
determining the number of D* and D° candi- 
dates in each VPDL bin. 

— The values of the parameters which had been 
fixed from the fit to "all" events, were varied 
by ±3<7. 

— The default bin width for the fits in the VPDL 
bin is 0.020 GeV. We lowered the bin width 
to 0.016 GeV and increased the bin width to 
0.027 GeV, and included the resulting varia- 
tions in the systematic uncertainty. 



X. CONCLUSIONS 

We have performed a study of a likelihood-based 
opposite-side tagging algorithm in B° and B + samples 
obtained with ~ 1 fb _1 of Runll Data. The dilutions 
V{B + ) and V(B a ) were found to be the same within their 
statistical uncertainties. This result justifies the applica- 
tion of the B® dilution to the B® mixing analysis. 

Splitting the sample into bins according to the tagging 
variable \d\ and measuring the tagging power as the sum 
of the individual tagging powers of all bins, we obtained 
a tagging power of 

eV 2 = [2.48 ±0.21 (stat.^oioe 0y st )] %■ 

From the simultaneous fit to events in all \d\ bins we 
measured the mixing parameter: 

Am d = 0.506 ± 0.020 (stat) ± 0.016 (syst) ps~\ 

which is in good agreement with the world average value 
of Am d = 0.502 ± 0.007 ps" 1 0. 
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TABLE IV: Systematic uncertainties for Am^. 





Default 


Variation 


Am d (ps x ) 






(a) 




(b) 


(a) 




(b) 


Br(B° -> D*~ii + v) 


5.44 


-0.23 




0.23 


0.002 




-0.002 


Br(B -> D'-k^vX) 


1.07 


-0.17 




0.17 


-0.0078 




0.0078 


R" 


0.35 







1.0 


0.0006 




-0.0012 


B lifetimes 


0.05022 


-0.00054 




0.00054 


0.0008 




-0.0008 


Resolution scale factor 




1.2 




0.8 


0.0021 




-0.0021 


Alignment 




— 10 /xm 




+10 ^m 


-0.004 




+0.004 


if-factor 




-2% 




+2% 


0.0098 




-0.0094 


Efficiency 




-12% 




+ 12% 


-0.0054 




0.0052 


Fraction D° in D* 


4% 


3.15% 




4.85% 


-0.0020 




+0.0030 


Fit procedure 


See below 


Bin width 


2 MeV 


1.6 




2.67 


0.0009 




0.0014 


Parameter fio 




-3a 




3a 


-0.0001 




0.0001 


Parameter aR '^' 7L 




-3a 




3a 


-0.0001 






Parameter <7R ~' T - L 




-3a 




3a 


-0.0001 




0.0001 


Parameter /ii 




-3a 




3a 


-0.0016 




0.0015 


Parameter ^±2^ 




-3a 




3a 


-0.0006 




0.0006 


Parameter Ft 




-3a 




3a 


-.0005 




0.0004 


Parameter (/12 — /J.i) 




-3a 




3a 


0.0006 




-0.0007 


Parameter 2J^£2, 




-3a 




3a 








Fit procedure 




Overall 


+0.0023 
-0.0019 


Total 






±0.0158 
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TABLE V: Systematic uncertainties for T>(B°). 











T>(B°) 


V(B a ) 


V(B a ) 


V(B a ) 




V{B a ) 






Default 


Variation 


0.1 < \d 


< 0.2 


0.2 < \d\ < 0.3 


0.3 < \d 


< 0.45 


0.45 < \d\ < 0.6 


0.6 < \d\ < 1.0 






(a) 


(b) 


(a) 


(b) 


(a) 


(b) 


(a) 


(b) 


(a) 


(b) 


(a) 


(b) 


Br(B° -> D*~ij + v) 


5.44 


-0.23 


0.23 








-0.001 


0.001 




0.001 


-0.001 


0.001 


-0.001 


Br(B -> D*ir^vX) 


1.07 


-0.17 


0.17 


0.0004 


-0.0004 


-0.0011 


0.0011 


-0.0019 


0.0021 


-0.0020 


0.0021 


-0.0008 


0.0028 




0.35 


0.0 


1.0 


-0.0009 


0.0016 


-0.0027 


0.0048 


-0.0042 


0.0079 


-0.0057 


0.0105 


-0.0066 


0.0124 


S lifetimes 


0.05022 


-0.00054 


0.00054 




-0.0001 


0.0001 


-0.0002 


0.0003 


-0.0001 


0.0003 


-0.0003 


0.0014 


-0.0003 


Resolution function 




xl. 2 


xO.8 


0.0005 


-0.0006 


0.0010 


-0.0012 


0.0020 


-0.0021 


0.0024 


-0.0028 


0.0028 


-0.0032 


Alignment 




— 10 [im 


10 (im 


-0.004 


0.004 


-0.004 


0.004 


-0.004 


0.004 


-0.004 


0.004 


-0.004 


0.004 


iC-Factor 




-2% 


+2% 






-0.0001 






0.0001 


-0.0001 








Efficiency 




-12% 


+ 12% 


0.0006 


-0.0007 


-0.0008 


0.0006 


-0.0012 


0.0011 


-0.0013 


0.0010 


-0.0021 


0.0019 


Fraction D° in D* 


4% 


3.15% 


4.85% 




0.0010 


-0.0010 




-0.0010 


0.0010 


-0.0010 


0.0010 


-0.0010 


0.0010 


Fit procedure 


See split below 


Bin width 


2 MeV 


1.6 


2.67 


-0.0026 


0.0002 


-0.0024 


0.0014 


-0.0001 


0.0027 


0.0037 


0.0038 


0.0089 


0.0087 


Parameter (j,q 




-3a 


3a 


-0.0003 


0.0002 


0.0001 


-0.0001 


0.0001 


0.0001 


-0.0002 


0.0001 


-0.0007 


0.0007 


Parameter €TR ^ €TL 




-3(7 


3a 


0.0002 


-0.0002 


0.0001 


-0.0001 


0.0004 


-0.0003 




-0.0001 


-0.0002 


0.0001 


Parameter a R , - L 




-3ct 


3a 


-0.0005 


0.0005 


0.0002 


-0.0001 


0.0002 


0.0001 


-0.0002 


0.0001 


-0.0015 


0.0011 


Parameter (X\ 




-3ct 


3a 


-0.0009 


0.0010 


-0.0017 


0.0018 


0.0023 


-0.0015 


0.0006 


-0.0005 


-0.0004 


-0.0004 


Parameter CT1 2 CT2 




-3ct 


3a 


0.0008 


-0.0005 


0.0014 


-0.0009 


0.0037 


-0.0034 


-0.0013 


0.0017 


-0.0099 


0.0068 


Parameter i? 




-3<t 


3a 


0.0015 


-0.0011 


0.0029 


-0.0024 


0.0030 


-0.0027 


0.0013 


-0.0011 


-0.0046 


0.0035 


Parameter (/Lt2 — Ml) 




-3a 


3a 




-0.0003 


0.0008 


-0.0011 


-0.0001 


0.0006 


-0.0003 


0.0002 


0.0008 


-0.0003 


Parameter ai . a2 
CT i 




-3a 


3a 


-0.0001 




-0.0004 


0.0003 


0.0002 


-0.0002 


-0.0004 


0.0004 


-0.0006 


0.0010 


Fit procedure 




Overall 


+0.0021 


+0.0040 


+0.0060 


+0.0044 


+0.0119 










-0.0031 


— 


3.0041 


-0.0046 


-0.0019 


-0.0111 


Total 








+.0049 


+ .0077 


+.0111 


+ .0125 


+.0182 










-0.0052 


— 


3.0066 


-0.0081 


-0.0081 


-0.0140 



